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ABSTRACT 

We determine the maximum lifetime t„ax of 52 FRII radio sources found in 26 central group galaxies from 
cross correlation of the Berlind SDSS group catalog with the VLA FIRST survey. Mock catalogs of FRII 
sources were produced to match the selection criteria of FIRST and the redshift distribution of our parent 
sample, while an analytical model was used to calculate source sizes and luminosities. The maximum lifetime 
of FRII sources was then determined via a comparison of the observed and model projected length distributions. 
We estimate the average FRII lifetime is L5 x 10^ years and the duty cycle is ^ 8 x 10^ years. Degeneracies 
between t„„u and the model parameters: jet power distribution, axial ratio, energy injection index, and ambient 
density introduce at most a factor of two uncertainty in our lifetime estimate. In addition, we calculate the 
radio active galactic nuclei (AGN) fraction in central group galaxies as a function of several group and host 
galaxy properties. The lifetime of radio sources recorded here is consistent with the quasar lifetime, even 
though these FRIIs have substantially sub-Eddington accretion. These results suggest a fiducial time frame for 
energy injection from AGN in feedback models. If the morphology of a given extended radio source is set 
by large-scale environment, while the lifetime is determined by the details of the accretion physics, this FRII 
lifetime is relevant for all extended radio sources. 

Subject headings: cooling flows — galaxies: active — galaxies: clusters: general — galaxies: evolution — 
galaxies: jets — radio lines: galaxies 



1. INTRODUCTION 

The interaction between active galactic nuclei (AGN) and 
their environment has been increasingly studied as a plausi- 
ble explanation of several long standing problems in our un- 
derstanding of galaxy formation and evolution. Hierarcha l 
structure formation, first formalized bvl White & ReesI (Il978l) , 
predicts a halo mass function that does not match the shape 
of the observed galaxy luminosity function (LF). Simula- 
tions predict an overabundance of massive and dwarf galaxies 
when compare d to the ob servable universe (IWhite & Fren3 
1199 It [Benson et alJl2003l) . Another conflict between obser- 
vation and theory is that galaxies are primarily red, early 
types with little star formation or actively star forming blue, 
late-ty pes (e.g. Blanton et al. 2003; Kauffmann et al. 2003, 
|2004|) . This bimodality requires abrupt truncation of star for- 
mation. Most conspicuously, the temperature distribution of 
the intergalactic medium (IGM) of so me clusters is at odds 
with si mple radiative cooling arguments (iPeterson et al.l200ll : 
iHicks & Canizares 2001). The explanation of all of these 
problems favors regulation of galaxy evolution by an un- 
known mechanism of AGN 'feedback'. 

The strongest case for the existence of AGN feedback 
can be made in the cores of galaxy clusters and groups. 
Clusters have long been known as str ong X-ray sources 
iByram et al. 1966; iGursky et alJ 11971 ). emitting thermal 
bremsstrahlung r adiation from their hot gas (first propo sed 
and confirmed by lFelten et"anil966l;lMitcheU et al.lll976l re- 
spectively). Calculations of the gas cooling time in s ome clus- 
ters showed tha t it was less than th e Hubble time (|Lea et alj 
[1973; Silk 1976). Subsequently, Cowie & BinneV ('1977[) 
and Fabian & Nulsen. (1977,) developed steady-state "cooling 
flow" models in which the gas of these systems condenses to- 



wards the center of the potential well. Two major problems 
with the cooling flow scenario have surfaced: (1) There have 
been no observations of the large star formation rate neces- 
sary for the center of clusters t o be mass sinks of the con- 
densing gas (Fabian et alJll99l1) and (2) High resolution X- 
ray observations confirm tha t cluster centers have far less gas 
below 1 keV than predic ted (Peterson et al. 2 0011; [Allen et al] 
1200 U iHicks et alJ [20021: [Fabian et al.. .2003.) . As evidence 
mounts that gas in these systems do not exhibit strong cooling 
flows, they h ave recently been renamed "cool core" clusters 
(iMolendi & Pizzolato 2001; Donahue & Voit 2004). 

While cooling flows have been discounted, galaxies ex- 
ist, so baryons must cool. Removing radiative cooling in 
the intracluster medium (ICM) produces the wrong s lope of 
the X -ray luminosity-temperature relati on (Ly - Ty ) (iKaised 
Il986t [Borgani et al. 2001; Muanwong eTal] l2001l) . To re- 
solve these discrepancies, radiative cooling must be balanced, 
at least in part, by some heating source. Initially, feed- 
back from supernovae was thought to heat the ICM, however, 
it is now assumed that only AGN could produc e the enor- 
mous energy needed to combat cooling (e.g. Vala geas & SiUd 
Il999t IWu et"an [20Q0). Early evidence for an AGN connec- 
tion camejrom cross correlating clu ster samples with radio 
sources (IOwen|[ 19741; lBahcal]|[l974h . Subsequent study by 
[lumi([i990F found that more than 70% of cluster dominant 
galaxies in cool core clusters are luminous radio galaxies. 
This radio-X-ray correlation led to growing interest in the pos- 
sibility of AGN feedback manifested as ICM heating via ra- 
dio lobes (e.g. [B innev & Tabor 1995; Chura zov et aj J 120011: 
ISoker et al.ll200k.Bruggen & Kaiser..2002;.Begelmanll2004 . 

Spectacular support of AGN feedback on the ICM came 
from the observations of large X-ray cavities in the gas sur- 



2 



Bird, Martini, & Kaiser 



rounding group and cl uster cores spatially co i ncident with 
the ra dio lobes of jets (iMcNamara et aljQoOOt iFabian et"an 
I2OOOI) . The prototype of radio mode feedback in clusters 
is now the set of extremely deep Chandra observations of 
Perseus, which show large cavities in the gas near the clus- 
ter core and e vidence for shocks , ripple s, and sound waves 
in the ICM (iFabian et all I2003L l2006h . Observations of 
Cygnus A also exhibit clear interactions between the cen- 
tral radio source and th e hot gas in this group of galaxies 
(e.g.lMionitir'2000^. ICM heating via radio mode feed- 
back would eliminate the need for a centr al mass sink and 
create the "entropy floors" seen in clusters (IValageas & Silkl 
Il999t iNath & Rovchowdhurv"2002V producing the observed 
Lx-Tx slope (Evrard & Henry 1991; Kaiser 1991). Most con- 
vincingly, the amount of energy necessary to heat the ICM 
is consistent with the energy needed to physically create the 
cavities in many systems (Birzan et al. 2004). AGN feedback 
probably occurs over a broad range of energies and efficien- 
cies: the sizes of a few radio cavities suggest that their cre- 
ation isjiot sufficient to completely balance radiative cool- 
ing dPope et al. 2007). The actual heating mechanism em- 
ployed by radio lobes is still under intense in vestigation (e.g. 
iRuszkowski et al.ll2004t [Reynolds et al.ll2005h . 

The timescale of AGN feedback is unknown yet fundamen- 
tal to understanding its underlying physics and to model ICM 
evolution and the accretion history of AGN. Several differ- 
ent methods have been introduced to estimate the AGN feed- 
back lifetime via th e radio emission of lobe s or X-ray obser- 
vations of the ICM. IOmma & Binneyl (l2004l) assumed that ra- 
dio cavities move at the sound speed of the ambient medium 
and calculate their age using their projected distance from the 
galaxy core. The time necessary for cavities to buoyantly rise 
to their prese nt position can provid e an estimate of their life- 
tim e as wen (IChurazov et al.ll2001h . iMcNamara et all (l2000l) 
andl Nulsen et al.l ( 2002h determined cavity age using the time 



required to refill the volume of ICM displaced by the rising 
bubbles. Cavity ages between ^10^ an d ~ 10^ years hav e 
been recorded using these techniques (e.g. lBirzan et al.l2004l) . 
However, some of these age estimates apply to the current 
AGN age, or the time since the AGN turned off. The analy- 
sis herein estimates the average total lifetime of FRII sources 
and consequently the timescale of their AGN energy injec- 
tion. This energetic AGN lifetime is critical to determining 
the physical processes of AGN feedback and whether radia- 
tive cooling can be balanced by said mechanisms. 

Radio loud AGN are not in every cool core cluster; 
thus, if radio-mode f eedback is viable, it must be episodic. 
iBinney &Tabo3(fT995.) were among the first to describe evo- 
lutionary models of the ICM. Radiative loses and gravitational 
forces promote gas infall towards the potential well and pro- 
vide fuel for the AGN. The AGN then forms jets which in 
turn lead to the creation of lobes that heat the ICM and cause 
a net outflow of gas. Once this fuel supply is removed, the 
AGN can no longer energetically sustain the jet and the newly 
created cavities then rise buoyantly through the ICM. The 
strongest evidence for episodic heating has been the discovery 
of "ghost" cavities in the outer regions of some clusters coin- 



ated with the AGN ( 


McNamara et al.1 120011 iJohnstone et alj 


120021 iMazzottaet all 


20021). 



feedback provides a natural regulation mechanism: the energy 
injected into the interstellar medium (ISM) prevents large 
scale star formation and provides an upper limit to the stellar 
mass of the galaxy, contributing to galaxy downsizing (e.g. 
ICowie et anil996t IScannapieco et a l. 2005). To better un- 
derstand the consequences of AGN feedback, hydrodynamic 
models have begun detailing the heating a ffect of radio lobes 
in the group and cluster environment (e.g. lSiiacki & Springell 
I2OO6I) . Hierarchal structure formation models now routinely 
invoke AGN feedback to obtain the observed shape of the 
galaxy luminosity function and bimodal galax y distribution 
e.g. Granato et al. 2004; Scringel et al. 2005aR [Bower et alJ 
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The preponderance of evidence supporting a connection be- 
tween galaxy evolution and AGN has led to the assumption 
that black hole growth regulates star formation ( Silk & Rees 
[19981 iFabiadll 9991: lKingi2(M I Weinmann et al.ll2006l) . AGN 



2006i: iCroton et aL.200fe .Hopkins et aL.2006) . In addition to 
the nature of the physical process that couples the radio source 
to the ICM, the two main sources of uncertainty in the afore- 
mentioned models are the efficiency of energy injection and 
again, the timescales over which the process occurs. We focus 
on the latter unknown to provide context for the consequences 
of AGN feedback and its role in galaxy evolution. 

Some new insights into the details of the dynamics of ra- 
dio lobes and how they transfer energy to the surrounding 
medium can be gained from high-resolution observations of 
individual sources. Multicomponent radio sources are usu- 
ally divided into one of two morphological classes : Fa- 
naroff and Riley Class I and II jFanaroff & Rilev[[T974b . FRI 
sources are dominated by emission towards the nucleus with 
a gradient towards fainter and more diffuse emission as one 
moves radially outward. FRJIs are limb brightened sources, 
with well defined areas of brightest emission ("hotspots") 
in their outermost regions. In addition to their morphologi- 
cal dis tinction, FRIIs ar e generally more luminous than FRIs 
( jFanaroff & Rilevlll974l) while the most powerful FRIIs usu- 
ally show a strong A GN cornponent in their spectra indicative 
of obscured quasars (lBarfliellll989l:[Hes et alJl993i) . FR class 
may also be a function of environment as the majority of ex- 
tende d radio sources in clusters are FRIs (Prestage & Pe acockl 
[19881: iHiU & LiUyllT99ll) . FRI sources may be FRIIs that en- 
countered high density environments as they arise from turbu- 
lent disruption of their lobes and entrainment of ambient gas. 
Two morphological groups within the FRI distinction. Nar- 
row Angle Tails (NATs) and Wide Angle Tails (WATs), are 
the result of radio source interaction with the ICM in dense 
cluster and group environments, either disp laced du e to galac- 
tic mo tion relative to the ambient m edium ([Owen & Rudnic^ 
[19761) or bulk motions in the ICM ([Burns et al.[[1994t) . Bofli 
FRI and FRII sources have been observed intera c ting with the 
surrounding mediumje.g. [Blanton et alJ [20001 : [Fabian et alJ 
,2003; .Wilson et al...2003l) . 

The weU defined structure of FRIIs is an astrophysical 
analog of a supersonic, collimated flow into ambient gas 
([de Young[ 12002). Long identified as the collision site of 
the energetic jet particles and the surrounding medium, the 
hotspot position of FRII sources has been used as a stan- 
dard measuring stick of the spatial extent of these powerful 
radio sources and can potentially reveal source age. Asymme- 
tries in the length of classic double lobed sources due to light 
travel time have been u s ed to e stimate their expansion speed 
([Longair & Riley[ [1979[: [Arshakian & Lon gair 2004). High 
resolution measurements of these geometric asymmetries es- 
timate the hot spo t advance spee d to be ~ 0.02c throughout 
the jet's lifetime ([Scheueil[1995h . The radiative properties 
of radio lobes can provide estimates of the hotspot advance 
speed and source age as well. Radiating particles advect away 
from the hotspot and fill the cocoon, losing energy and soft- 
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ening the radio spectrum of the source as a function of age 
(lAlexandeiill987HAlexander & Leahylll987h . Detailed obser- 
vations of Cygnus A show that the spectral inde x is a stron g 
function of radial distance from the AGN (Carilli et al.lll99lT) : 
this is seen in other FRII sources as well. Recent spectral 
age estimates of FRIIs are between 10^ and 10^ years (O'Dea 
2007, in prep.). 

Shared characteristics of FRII sources led to the idea of self- 
similar growth (Falle 1991). In addition to similar advance 
speeds, similar axial ratios (Rt, the ratio of a lobe's length 
to width) were recorded for F RII sources over a wide range 
of linear sizes ( Leahv & Williimi Il984t iLeahy et all 119891: 
ISubrahmanyan et al.al996i) . According to the basic model for 
FRII growth, an A GN gives rise to two jets that propagate in 
opposite directions ('Blan dford & Reeslll974l: IScheuer" 1974j). 
The energetic jet terminates in a shock, forming the hotspot, 
and drives a bow shock into the ambient medium. This ac- 
cepted foundation has been built on by a number of mod- 
els th at analytically describe the dynamics of FRII sources 
(e.g. Kaiser & Alexander 1997, hereafter KA) and their lumi- 
nosity evolution (e.g. Ka iser et al. 1997; Blundell et al. 1999; 
iManolakou & Kirkll2002i hereafter KDA, BRW, MK, respec- 
tively). All three use the dimensional arguments of KA to 
determine the length of FRII sources as a function of age but 
differ on how energy is deposited in the radio lobes from the 
hotspot. Given the physical properties of the jet and its ambi- 
ent environment, each model predicts the track of sources in 
the luminosity-size plane. We use the KDA model as the cru- 
cial link between the intrinsic and observable characteristics 
of FRIIs. While the physics of radio jets are not completely 
understood in detail, these models of FRII source evolution 
produce res ults consistent with th e demographics of large ra- 
dio surveys (iBarai & Wiitall20()6h . 

Demographic studies of radio mode feedback have been 
greatly assisted by the recent completion of large area op- 
tical and radio sky surveys. In particular, the Sloan Dig- 
ital Sky Survey (SDSS), covering two fifths of the North- 
ern Galactic sky down to an r-band limiting magnitude of 
^ 22.5, has proved invaluable in determining where poten- 
tial radio mode feedback systems (groups and clusters) are on 
the sky. Several systematic searches for clusters and groups 
using photometric (e.g. lAnnis et al.l II999I ; iKim et al.n2002l; 
Goto et alJ2002HLee et alJ2004 and spectroscopic (e.glGotol 
2005t iMiUer et aljT2005l: iMerchan & Zandivare3l2005h data 
have provided valuable catalogs for a plethora of research. 
The Faint Im ages of the Radio Sky at Twenty centimeters 
(iBecker et alJll995l: IWhite et al.lll997i FIRST) survey is the 
radio analog to SDSS. Surveying the SDSS area down to 1 
mJy at 1 .4 GHz, FIRST represents a substantial improvement 
in radio surveys with regards to sensitivity and angular reso- 
lution. Ot her radio surveys, especially the NRAO VLA Sky 
Survey dCondon et a l. 1998, NVSS), and FIRST constitute 
an expansive and detailed picture of the radio universe. The 
combination of these optical and radio resources provides un- 
paralleled opportunity to investigate the radio properties of 
the richest systems in the universe. In addition to our FRII 
lifetime estimate, we present a demographic study of radio 
sources in low redshift central group galaxies. 

One of the first investigations from these large area sur- 
veys has been the measurement of the ra dio galaxy LF and 
its deviations fro m the optical galaxy LF (ISadler et al.ll2002l : 
iMauch & S adler 2007). Correlations between radio active 
AGN and host galaxy properties, including local density and 
luminosity, have strengthened the arguments for radio mode 



feedback (i Best et aT]|2005allbl) . In addition, a higher fraction 
of radio loud AGN are found in galaxies at the center of clus- 
ters compared to simila r galaxies not at the center o f their 
cluster's potential well dBest et al.l [20061: ICroft et al.l 12007). 
Demographic studies have provided good evidence that ra- 
dio mode feedback is a factor in galaxy evolution and even 
regulation. 

In the present study we produce a sample of FRII radio 
sources associated with central group galaxies using optical 
data from the SDSS and radio maps from FIRST to measure 
the maximum lifetime of FRIIs. Our systematic search for and 
subsequent analysis of FRIIs measures their typical energetic 
lifetime- a fundamental parameter necessary to ascertain the 
nature of AGN feedback and its ramifications. The properties 
of the host galaxy are assembled from our optical data while 
the source attributes are obtained from our radio data. We ul- 
timately produce a collection of FRII radio sources and mea- 
sure their projected lengths (Section|2]). Simulations of FRII 
sources with known properties are produced in a Monte Carlo 
scheme using the KDA model and compared to our observed 
sample in Section [3] In Section |4] we present our results for 
the lifetime of FRII sources, while dependence upon group 
and host galaxy properties is discussed in Section|5] We sum- 
marize our findings in Section|6l 

2. SAMPLE 

2.1. Optical Sample 

This study's primary sa mple is a subset of the group and 
cluster catalog created by iBerlind et alJ (l2006ah using the 
SDSS. Using a dedicated 2.5 meter telescope at Apache 
Point Observatory, SDSS scans the sky in five photomet- 
ric bandpasses simultaneously with its 120 Megapixel cam- 
era dGunn et al.l 11998b . Except in the case of fiber colli- 
sions owing to very small angular separations on the sky, 
most ga laxies in the survey ar e selected for spectroscopic fol- 
low up. iBerlind et al.l (l2006ah use a redshift space friends of 
friends (FoF) algorithm to identify groups of g alaxies within 
the SAMPLE 14 galaxy set (iBlanton et al.ll2003l); itself a sub - 
set of the main SDSS galaxy sample (IStrauss et ani2002h . 
FoF identification schemes are based on the premise of re- 
cursively linking galaxies together within a specified link- 
ing volume around each galaxy. The specific details can 
be found in Ber Und et al.l (2006a) and sources cited therein. 
FoF algorithms have several advantages: they produce unique 
groups of galaxies for a given linking volume, no group ge- 
ometry (e.g. spherical) need be assumed a priori, and groups 
are only supplemented with new members when the linking 
volume is increased. The Ber lind et al.l (2006a) group cata- 
log is a volume limited sample of 57138 galaxies complete 
down to Mr = -19.9 mag and spanning the redshift range 
0.015 < z < 0.1. Each group member's luminosity and spec- 
troscopic redshift, in addition to the total group luminosity, 
are reported in the catalog. 

2.2. Radio Sample 

We used the FIRST (iBeckeret al.llI995h survey conducted 
at the VLA to identify radio sources in our group sample. 
FIRST, a radio survey conducted at 1 .4 GHz, has an angular 
resolution of ^ 5" and a typical RMS of 0.15 mJy per pixel. 
The resulting catalogs have a source detection limit of 1 mJy 
and astrometry better than 1". The survey's area was designed 
to match that of SDSS. FIRST represents roughly an order of 
magnitude improvement in angular resolution and sensitivity 
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over previous large area radio sky maps. These survey prop- 
erties allow us to identify extended radio sources associated 
with single galaxies even in a crowded group environment. 

FIRST observations of FRJI sources were complimented 
with data from the NRAO VLA Sky Survey (NVSS). Also 
carried out at 1.4 GHz, NVSS is complete do wn to 2.5 mJy 
and h as an angular resolution of 45 arcseconds dCondon et alj 
Il998h . Due to NVSS's shorter baselines in comparison to 
FIRST, it is much less likely to resolve out any flux from 
sources. 

2.3. Cross-Matching of the Optical and Radio Catalogs 

To create a manageable sample size, groups with total 
Mr < -22.0 mag were included in this study. In addition, 
groups whose central galaxy's coordinates were not included 
in FIRST were excluded. The remaining sample is composed 
of 2020 groups with mass greater than ^ 10'^'^/i~'Mq and 
a typical velocity dispersion range of 100 < ^ .700 km 
s"' . Central group galaxies, identified in lBerlind et alj ( l2006ah 
as the brightest group member, were cross-correlated with 
FIRST. 

We note that the FIRST catal og reports radio "co mponents" 
as opposed to whole sources (IWhite et alJII997)). As such, 
any FRII source would have two components (Ide Vries et alJ 
I2OO6I) . We searched for multicomponent (> 2) radio sources 
in the FIRST catalog within 100 kpc of each central galaxy. 
Each candidate was then visually inspected for FRJI morpho- 
logical features (limb brightened double sources). FR mor- 
phology can be robustly detected when the source is > 10" 
(corresponding to the minimum FIRST component size plus 
the approximate half power beam width of FIRST) and we re- 
quire this minimum length of any potential FRII. We find 32 
central galaxies whose associated radio sources meet this sim- 
ple criteria. Any FRJ interlopers are removed from the sample 
via the morphological distinction made by Fanaroff & RilevI 
(|T974), i.e., if the "ratio of the distance between the regions 
of highest brightness of opposite sides of the central galaxy 
or quasar to the total ext ent of the source" was g reater than 
0.5, the source is an FRII dFanaroff & Rilevlll974 . Of our 32 
potential candidates, 26 extended radio sources pass this mor- 
phological cut. These galaxies and their 52 associated lobes 
comprise our sample and their properties are summarized in 
Table[T] Thumbnails of the optical and radio emission of each 
galaxy in our sample are shown in Figure [T] We discuss de- 
termination of lobe and "hotspot" length in Section 12.41 

After the initial FR morphological distinction by 

lFanaroff& RilevI ( fl974l) . FR sources have historically 
been classified using an empirically seen brea k between 
the luminosity of the tw o classes ( Fanaroff & RilevI 119741 : 
iLedlow&Oweiil 119961) . iLedlow & Owem 0996) show the 
break is a function optical luminosity {Lrad oc Lgpt) but is less 
well-defined in clusters and at lower redshifts. The mean M,- 
of the host galaxies in our sample is -21.82, corresponding to 
an FRI/FRII break luminosity of Liachz ~ 4.7 x lO^"* W Hz"' 
at z = 0.1 for our choice of cosmology (20 sources do not 
exceed this threshold in our sample, of which 9 are within a 
factor of ~ 3). We note that studies of FR sources have been 
dominated by high luminosity sources due to the relative 
lack of sensitivity in wide-a rea radio surveys until recent 
years. iBlundell et all (IT999I) argue that morphology is a 
much more fundamental quantity than luminosity as similarly 
luminous sources will not have the same jet properties at 
all redshifts. Our sample is strictly composed of FRJIs in 
the morphological sense. Further, KA and KDA model the 



self similar evolution of limb brightened lobes with distinct 
working surfaces. Our sample meets this criteria and thus the 
KA and KDA models can be applied appropriately. 

2.4. Projected Source Length Measurement 

Several methods have been employed in the literature to 
measure the projected length of radio lobes. In the past, a few 
studies investigated the ratio of projected lobe lengths of clas- 
sic double lobed radio sources to determine their e xpansion 
velocities (iLongair & RilevI 119791; [Scheuedll995l) . IScheued 
(1 1995!) measured the projected lobe length as the distance be- 
tween the "second or third" radio contour at the edge of the 
lobe and the core radio source, noting "that any well-defined 
prescription, rigorously applied, will sometimes produce ob- 
vious nonsense." In more recent work, the length of the source 
is defined as the distance from the "hotspot," or the termina- 
tion shock of the jet, to the central radio source. 

In order to securely define a sample of FRII sources, we 
must know both the lobe length and the hotspot to central 
galaxy separation. We assumed that the centroid of the bright, 
compact emission in each lobe was the position of the hotspot 
(as in Blundell et al. (1999)). This definition is practical as 
it traces the position of the brightest emission in each lobe 
and provides half of the needed inform ation to make FR class 
distinctions of iFanaroff & Rileyl (Il974l) . Determination of the 
lobe length is inherently more subjective. We define the lobe 
length as the maximum separation between the central galaxy 
and the contour corresponding to 5 times the average RMS 
of the FIRST survey (5 *0.15 mJy bm~'). The errors asso- 
ciated with both the position of the brightest emission and 
the lobe length defined in the above manner are less than one 
FIRST pixel (1.8"). Lobe lengths are measured exclusively 
using FIRST data as our study is dependent upon accurate 
astrometry to establish galaxy-lobe associations and to make 
source length determination as rigorous as possible. The cu- 
mulative projected length distribution of the observed sample 
has a median of 36.5 kpc and is shown in Figure |2l See Ta- 
ble[T]for the projected lobe lengths of the FRII sample. 

In our subsequent analysis we only include lobes with a 
projected length less than 100 kpc. This cut in lobe length 
removes five lobes from our FRII sample but suppresses two 
significant sources of potential error: (1) Longer sources in- 
herently have larger absolute errors in their length measure- 
ment and (2) At distances greater than 100 kpc from the cen- 
tral group galaxy, the potential for associating another group 
member's radio emission with the central galaxy increases. In 
addition, we expect at least one convincing fake FRII source 

o 

(greater than 135 separation between two random sources 
with projected length asymmetry less than 15 kpc) at greater 
distances due to the surface density of FIRST sources. We 
examine any effect this cut may have on our conclusions in 
section l4l1 

3. MOCK POPULATION GENERATION 

To estimate the maximum age of the sources in our sample, 
we generate mock FRII radio source catalogs that match the 
selection criteria of our FIRST sample but span a wide range 
of maximum age. We then compare the model distributions 
of lobe length with our observed distribution to determine the 
maximum lifetime imax most consistent with the data. There 
are several considerations one must take into account when 
creating the mock catalogs: the random projections of our ob- 
served jets: the limitations of FIRST (both in angular reso- 
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TABLE 1 
FRII Source Sample 
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Note. — FRII sources in our sample inorder of increasing right ascension. (1): SDSS identifier. (2); Galaxy redshifi. (3); Total r-band absolute magnitude of the group. (4); (g-r) color of host group. (5); Absolute r-band magnitude of the host galaxy. (6); 
(g-r) color of the host galaxy. (7): Total source fiux in mJy reported by FIRST. (8); Total source flux in mJy reported by NVSS. (9); Total source luminosity in lO'^^ W Hz~^ computed from the NVSS flux. (10,1 1): The projected separation of the hotspot 
and central galaxy in kpc. (12,13): The angular size of the two lobes in arc seconds. (14,15); The projected length of the two lobes in kpc. Note that (10,12,14) will always correspond to the east hotspot and lobe, (11,13,15) to the west hotspot and lobe. 

lution and sensitivity); and the evolution in luminosity and 
linear size of radio sources as a function of time. 

A simple Monte Carlo scheme is employed to create the 
mock catalogs. For a given simulation, sources are assigned 
ages, redshifts, jet powers, and orientations according to the 
prescriptions in this section. We detail the age, redshift, 
and jet power distributions of our catalogs in Section 13.11 
Each source is then evolved using the self-similar hydrody- 
namic KDA model to produce their luminosities and intrinsic 
lengths. A brief outline of the KDA model and its input pa- 
rameters follows in Section 13.21 Sources are projected onto 
the sky with a random viewing angle (Section [3. 3l l. We as- 
sume a flat universe with Hq = 71 km s~', VLm = 0.27, and 
f7A = 0.73. 



3.1. Global Population Properties 

In order to model the luminosity and length of a mock ra- 
dio source, the three input parameters are the age, redshift, 
and jet power. We generate a series of mock catalogs parame- 
terized by the maximum allowed age l„jax of the sources. Each 
mock source is assigned some time, tj, between its birth time 
and l,„ax- For a single simulation's sources, t^- is uniformly 
distributed over the interval [1 yr, \.„,ax}- In practice it is rea- 
sonable to assume some distribution in t„„u as each source 
will not turn off at exactly the same time. While one could 
in principle attempt to solve for a spread in t,„nv as well, we 
do not do so because of our small sample size. Introducing a 
distribution in imax to the model would widen the likelihood 
functions calculated in Section H] Once a source is assigned 
tj, it is then randomly assigned a redshift from a probability 
distribution that matches the observed redshift distribution of 
the parent group sample. 



Finally, each source is assigned an initial jet power, Q, 
drawn from a probability distribution of the form: 

p{Q)dQ^Q-''dQ, (1) 

for Q in some range 2„„„ < Q< Qmax- We adopt two dis- 
tinct Q distributions, performing a parallel analysis with each. 
The first is const rained by the local LF of active radio galax- 
ies calculated bv lSadler et alj (120021) . Jet powers are assigned 
to our mock sources from a probability distribution such that 
when evolved with KDA under our default conditions (Ta- 
ble ID, the resultant luminosity function's slope and limits 
match that of radio emitting AGN (see eq. 11. ISadler et alj 
I2OO2I) . The Sadler Q distribution (hereafter Qs) designates: 
Qmin = 5.0 X 10^3 W, e„,„ = 1.0 X 10^ ^ W, and a = 62. Th e 
second Q distribution is taken from ' Blundell et all (119991) : 
Q,nin = 5.0 X 10" W, Q,„ax = 5.0 X lO'*^ W, and a = 2.6 (here- 
after Qb). While Qmax is higher for this distribution, the much 
steeper slope yields fewer sources with high Q. Neither of 
these Q distributions are ideal for our purposes: The Sadler 
Qs is derived relative to an FRI LF, while the Blundell Qb 
employed a different cosmology. However we demonstrate 
below is Section |4]and the Appendix that our results are rel- 
atively insensitive to the Q distribution. Note that sources 
within a single simulation are assigned jet powers from only 
one of these Q distributions. 

3.2. KDA Evolution Model 

The next ingredient in our mock catalogs is a prescription 
for radio source evolution in time. As noted above, three 
models for FRII source evolution are commonly discussed in 
the literature: KDA, BRW, and MK. Most of fliese models 
only differ in their calculation of the luminosity, whereas the 
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Fig. 1 . — The host galaxies with FRII radio sources in our sample. The SDSS r-band images are the grayscale while the contours ai'e logaiithmic radio flux 
from FIRST. Each image is 162" on a side. The scale bar in the bottom left is 30". 



equations governing the length of the radio source as a func- 
tion of age simply stem from a dimensional argument found 
in KA. iBarai & Wiital ( l2006l) quantitatively compared these 
models and their ability to reproduce established radio sur- 
veys' results, particularly the FRIIs in the Third Cambridge 
Revised Revised (3CRR), Sixth Cambridge (6CR), and Sev- 
enth Cambridge Redshift Surveys (7CRS). These authors cre- 
ated mock catalogs (using a similar scheme to this paper) to 
emulate these surveys and evaluated each model with ID and 
2D Kolmogorov-Smirnov tests between the real and modeled 
data across several parameters (luminosity, linear size, and 
spectral index). They concluded that while all of the mod- 
els required modification to produce the proper distribution 



of spectral index, the KDA model produced simulated data 
that most closely matched the observed surveys, especially in 
the luminosity-linear size plane most relevant to our investi- 
gation. Based on Barai & Wiita (2006), we use KDA to model 
the luminosity evolution of FRIIs in our mock catalogs. While 
we chose the KDA model based on its agreement with exist- 
ing data, all the models mentioned produce broadly consistent 
size-luminosity relationships and therefore our results should 
not significantly depend upon model choice. 

KA and KDA also have the virtue of being the simplest self 
similar models for the length and luminosity of a radio source 
as a function of age, environment, and jet power They are 
an extension of canonical models in which a jet emerges from 
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Fig. 1. — continued. 



the AGN region and is soon confined by the uniform pres- 
sure of the surrounding cocoon except near the hotspot. The 
ram pressure of the jet is distributed over the working surface 
and is balanced by the ambient shocked IGM. KA introduce 
the dynamics of the cocoon while KDA calculate evolution- 
ary tracks for radio sources in the luminosity-linear size plane. 
Below we quickly summarize the important assumptions and 
quantities involved in the calculation of source length. 

The KA model assumes the lobe expands into an IGM 
whose density is parameterized by a simplified King model: 
Px = Poir/aoT^ with r the radial distance from the radio 
source core, po a constant density, gq a scale length, and (3 
the radial density index (KA, see their eq. 2). Further assum- 
ing that the rate of energy injection from the AGN into the 
cocoon and rest-mass transport along the jet are constant, one 
can make a dimensional argument that the length of the lobe, 

is: 

.fx 3/(5-/3) 

Z = ciflo(^-j , (2) 



where 



Go 



(1/3) 



(3) 



is a characteristic time scale, t is the age of the source, Qo 
is the initial jet power, and ci is a dimensionless constant of 



order unity. The calculation of source luminosity is far more 
complex but still relies upon the basic assumptions outlined 
above. We use eq. 16 of KDA to compute the luminosity, P,y, 
as a function of time and frequency and refer the reader to that 
paper for its derivation. 

Relative to many radio source evolution models, the KDA 
approach is straightforward. Table|2]lists the parameters used 
in the model, a short description of each, and typical values. 
Given a choice of these parameters, hereafter collectively re- 
ferred to as and the age, redshift, and jet power of the 
source, the KA and KDA models describe the calculations 
necessary to obtain its luminosity and length. We first address 
the parameters that remained constant in our analysis and then 
discuss those that were varied from simulation to simulation. 

The adiabatic index of the gaseous IGM, cocoon, and mag- 
netic field (Fv, Fc, and F^, respectively) are chosen to repre- 
sent non-relativistic (Ft) or relativistic (F^, Tb) equations of 
state i.e. F^ = 5/3, F^ = Fb = 4/3. There are two other plau- 
sible scenarios: both the cocoon and the magnetic field are 
governed by non-relativistic equations of state or the magnetic 
energy density is proportional to one of the relativistic parti- 
cles. Neither scenario would significantly change the length 
calculation. 

Several jet parameters can float within reasonable ranges. 
The first of these is the energy distribution of injected parti- 
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TABLE 2 

Description and default values for KDA input parameters 



Fig. 2. — The cumulative distribution of projected lobe lengths in our 
sample (solid histogram) and representative models with tma_\= (dotted 
hne), lO' (solid line), and 10** (dashed line) years. The remaining parameters 
are our default values, X[s (see Table[3). These CDF shapes correspond to 
equally unique PDFs; thus, we are able to easily determine the most likely 
t;ira.v for our observed sample. Of all the variables in our analysis, t„,a.v has the 
largest effect on the length distribution. The best fit t„,o.v = 1.5 X lO' years. 

cles, which characterizes how energy is transported from the 
hotspot to the lobe. In the KDA model, radiating particles are 
injected into the lobe from the hotspot according to some en- 
ergy number distribution, N(E) cx E~p, and the injection index 
p is assumed to be constant over the lifetime of the source. In 
the rest frame of the jet shock, material in the jet is moving 
with r - 1.67 (KDA). Heavens & Drury (1988) determined 
that this implies an injection index of p = 2.14. While that ap- 
proach was not spec i fically intended for extragalactic shocks, 
[Alexander & Leahvl (Il987h and others have shown that typ- 
ical radio source SEDs imply 2 < p < 3. Therefore, KDA 
adopts p = 2 .14 as a fiducial value and we do the same. In 
section l4.2.3l we explore the relationship between p and t,„a.i, 
including the effect a more conventional index p = 2.5 has on 
the calculated maximum lifetime. 

The remaining parameters set the intrinsic properties of the 
lobe and the surrounding IGM. The axial ratio of the lobe Rj 
is the ratio of the length of the lobe to its width and is an 
important variable in the lobe's evolution both in power and 
as a linear source. The higher the axial ratio of the lobe, the 
higher the pressure at the hotspot (/?/,). This higher /?/, will 
produce a greater hotspot advance speed, yielding sources of 
a given length in a shorter period of time. The distribution of 
axial ratio in radio sources is still not well characterized by an 
unbiased su rvey of FRII source s. Work bv lLeahv & Williamsl 
(IT984i) and iLeahy et aT] (Il989h at least indicate axial ratios 
1.3 < Rt < 6. More recent samples confi rm this range as 
nominal for FRII sources (e.g.lKharb et al.'2007). Ideally, one 
would assign axial ratios to each model source from a known 
distribution. However, this distribution is poorly constrained 
by current observations and we therefore assign all sources a 
single axial ratio per simulation. KDA adopt a nominal value, 
Rt= 2.0, corresponding to = 31.1° and we do the same. 

A dense IGM increases the confinement pressure on the 
lobe and reduces the hotspot advance speed. We investigate 



Parameter 


Default Value 


Description 


r, 


5/3 


Adiabatic index of the IGM 


r. 


4/3 


Adiabatic index of the cocoon 




4/3 


Adiabatic index of the magnetic field energy density 


p 


2.14 


Injection index 


Rt 


2.0 


Axial Ratio 


PO 


7.2 X 10-22 kg m-3 


Constant central density 


ao 


2.0 kpc 


Scalelength 


13 


1.9 


Radial density power law index 



how the density profile of the IGM, another input parameter 
of the KDA model, affects lobe length and lifetime measure- 
ments in Section l4.2.1l and determine the relationship between 
axial ratio and t„,ax in Section l4.2.21 

3.3. Geometric Considerations 

As our actual measurement is of a distribution of projected 
source length, we finally assign each mock source a random 
viewing angle. This is more accurate than assuming an av- 
erage viewing angle for all sources. In the absence of any 
luminosity evolution, the expected source size distribution for 
a single age population, as well as for a range of ages up to 
some t,„aj, can be calculated analytically. For a single age 
population, where r is the intrinsic source length and R is the 
projected source length, the fraction of observed sources with 
length less than R is; 



f«R)=l-\ll-[^ 



(4) 



This function is shown in the left panel of Figure [3] along 
with the projected length distribution of 10** random sources 
of a single age. Each source is assigned a random d,(f) and 
projected on the sky, yielding the length distribution shown 
by the data points. Our random distribution is in excellent 
agreement with the analytic result. 

We now extend this test to a set of sources uniformly dis- 
tributed in age up to some maximum age t„ax associated with 
a maximum intrinsic length r,„aj . The fraction of sources less 
than some projected length R is: 



f«R)= / dr+ 



1- 




(5) 



We naturally split the equation between contributions from 
sources with intrinsic length r < R and r > R. For all r < R, 
/(< R) picks up the number of sources with intrinsic length 
r. For other values of r, we sum over the fraction of sources 
less than R for each intrinsic length. The right side of Fig- 
ure [3] shows /(< R) for a 10^ source sample (data points) 
overlaid with eq. |5] Notice how the projected length distri- 
bution is shifted towards smaller source lengths, with the me- 
dian source length now at 36% of r,„„x. In practice, this an- 
alytic distribution is not a good match to the data because of 
the implicit assumption of no luminosity evolution. Because 
sources fade as they age, the observed distribution should be 
skewed towards shorter sources. 

3.4. Catalog Generation 

For a given choice of parameters X (see Table |2]i and max- 
imum age t,„flj, we generate a mock catalog with lO*" radio 
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Fig. 3. — The cumulative fraction of sources less tlian a projected length R on the sky. The boxes are the result of random projections of single (left) and 
uniformly (right) aged sources. The maximum length is normalized to one in both plots. For a population of sources with a single age (left). 50% of the sources 
are projected at ~ 85% or higher of their intrinsic length. Once the population is uniformly distributed in age, the median projected distance is now 36% of the 
intrinsic length of the oldest sources. The boxes are the result of random projections of 10"* (left) and lO' (right) sources. The solid lines are given by eq.|4](left) 
and eq. [3] (right). 



sources. Each source is assigned a jet power, redshift, and 
age according to the distributions detailed above. The intrin- 
sic length and luminosity of each source is obtained using 
the KDA model. The se ar e randomly projected on the sky 
as described in section [331 yielding projected length Ip. The 
mock catalogs are stripped of undetectable sources via length 
and luminosity cuts to match the observational limits of the 
FIRST survey. If Ij, is less than our 10" minimum angular 
size for detected FRIIs, it is removed from the final catalog. 
Similarly, if P^, corresponds to a flux below the 1 mJy sensitiv- 
ity of FIRST, it is discarded. We examined a similar catalog 
generation process using a surface brightness cut and noted 
no statistical difference in the resultant mock catalogs when 
compared to those created using the FIRST flux limit. This 
demonstrates that our analysis is not biased against large, dif- 
fuse, presumably old lobes that might have been resolved out 
by first's exquisite angular resolution. The output from this 
process is a cleaned catalog of sources whose properties are 
well understood and meet our selection criteria. 

4. RESULTS: LIFETIME OF FRIIS 

If the physical properties of a set of radio lobes are kept 
constant, their maximum age sets the length distribution. In 
the next subsections we compare the projected length distri- 
bution of the models discussed in Section [3]with the observed 
distribution and quantify the goodness of fit via a maximum 
likelihood estimation (MLE). We first describe our maximum 
likelihood approach in Section 1431 and use this framework to 
estimate t„„u . We then investigate degeneracies between t^ax 
and other parameters of the KDA model in Section|42| 

4. 1. Method of Maximum Source Age Determination 

We use the maximum likelihood approach to determine 
which models best reproduce the length distribution observed. 
In general a model is made distinct by a unique choice of \.„u,x 
and model parameters (Section 13.21 and Table |2]l. Every 
choice of Xi and \.„u,x produces a distinct distribution of pro- 
jected lobe length I p. 

We first create probability density functions (PDFs) of the 



model distributions we wish to compare. For each mock cat- 
alog, we divide the one dimensional space of Ip into equally 
wide bins. The probability, Vk, assigned to the A:th bin is then 
just nk/J2j'^j where n^ is the number of sources within the 
^th bin and J^j^^i is the total number of sources in all the 
bins. The likelihood, C, of a model is then: 



(6) 



where Vj is the probability of finding the /'th source in the 
model. As the PDFs are relatively smooth functions, we use 
simple linear interpolation to obtain Vj from Vk- 

In order to reasonably limit the theoretical error of our 
model with a feasible number of sources, we only evaluated 
the probability of observed sources with 5.3 < I < 101 .3 kpc. 
Using these length restrictions, we are assured good statistics 
in our longest length bins and avoid potential sources of error 
such as spurious associations (Section [2.4b and sources that 
fall below the angular resolution of FIRST at the highest red- 
shift of the parent sample (z = 0.1). The maximum length of 
101.3 kpc was simply chosen to produce 33 equally spread 
bins of 3 kpc each. As we show below, these cuts in length 
do not limit our sensitivity to maximum age distinctions be- 
cause the observed length distribution is such an important 
constraint (see Figure O. We create 10* sources per simula- 
tion to ensure uniform sampling of the chosen distributions of 
jet power, age, and redshift. 

The prescription for determining the lifetime of FRII ra- 
dio sources in now straightforward. For a fixed choice of A/, 
we produce models with a range of t,„aj and determine that 
which maximizes the likelihood. Specifically, models with 
2 X lO'* sources were created for a range of t,„ax between 10* 
and 10*^ yrs in steps of 10* yrs. Once the approximate peak in 
the likelihood function is found, we produce models with 10* 
sources, making smaller steps in t,„a.x around the peak. The 
lifetime of FRII sources for a given A/ is then the t„„u of the 
model of maximum likelihood. One cr errors, corresponding 
to AlnC = 0.5, are calculated in each direction. Table[3]lists 
the lifetime of FRII sources for a variety of Xi. 
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TABLE 3 
Maximum Lifetime Estimates 



ID po ao P Rt p 1,„„(10'' yr) Notes 

(kgm-3) (kpc) Qs Qb 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 





7.2 X 10-22 


2.0 


1.9 


2.0 


2.14 


1 22±''28 

ii5±8:M 


1 75+0-21 
am 


default (KDA) 




L67 X 10-23 


10.0 


1.5 


2.0 


2.14 


1 59+'' " 


BRW density profile 




7.19 X 10-2<' 


391.0 


1.23 


2.0 


2.14 


1 31+''-"' 


Jetha et al. (20071 density orofile 




7.2 X 10-22 


2.0 


1.9 


1.0 


2.14 


^•^^=■=0.27 


5 ±8:i9 
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Figure |2] shows the cumulative length distribution of our 
sample, along with mock catalogs generated with our default 
KDA input parameters iX\s) for t„,a.,= 10^, 10^ , and 10^ years. 
It is important to note how dramatically different the model 
distributions are for different t,„ax- The fact that these distri- 
butions are so distinct allow us to determine t,„a_t with good 
precision and makes this technique so effective. For Xis and 
imax= 10^ yrs, the median projected lobe length is ~ 7 kpc, 
yielding no sources > 11 kpc. Conversely, a lifetime of 10*^ 
yrs produces of median source length of ^ 53 kpc. Shorter 
tmo.v skews the distribution towards young sources that do not 
have enough time to produce the longer observed lengths, 
whereas most sources are older than 10^ years in the t,„„v= 10** 
yrs case and therefore the expected observed distribution is 
much more heavily weighted towards longer lobes. Figure |4] 
shows the PDFs of several t,„ai for X[s. Even At„a^^ 2 x 10^ 
years produces appreciable changes in the PDF and thus the 
likelihood, C The t,„af best representing the data, the maxi- 
mum of the C{tmax) function, is obvious given a choice of X 

The results for our fiducial parameter sets, X\s.\b are shown 
in Tabled For X^s the most probable t„„^= 1.22(±jj;^^) x 10^ 
yrs, while for the Blundell distribution Xig the most likely 
tmfl.v= 1 -^SCibQ Q^) X 10^ yrs. A plot of the relative £ as a func- 
tion of \.„u,x for X\s and X\b is shown in the left panel of 
Figure |5] These two maximum age estimates are compara- 
ble because the median jet powers of the two Q distributions 
are similar in spite of their different ranges and slopes. The 
Blundell Q distribution has a slightly lower median jet power, 
hence the longer estimated lifetime of sources. 

4.2. Parameter Dependence 

We have determined the F Rll li fetime for a fiducial choice 
of model parameters (section 1431 1. An immediate concern is 
how degenerate the FRII lifetime is with various properties 
of the jet, lobe, and ambient IGM. We produce a series of 
models representing the extremes of the observed param- 
eter space discussed in Section [32] We estimate the lifetime 
of FRII sources with these parameter choices to place con- 
straints on the systematic uncertainties due to the unknown 
values and distributions of these parameters. As these mod- 
els should bracket the true range of jet, lobe, and environment 
properties, we expect the corresponding range of t„„j^ will be a 
conservative estimate of its uncertainties. See section l4!2.4l for 
a brief explanation of the choice of singular values over dis- 



tributions of Rt and density profile in a given simulation. An 
analytic investigation of parameter dependence is presented in 
the Appendix. 

4.2.1. Density Profile 

The density of the IGM surrounding a lobe could have a 
significant impact upon its length. This introduces a possible 
degeneracy between a chosen density profile of the IGM and 
the lifetime we measure for our FRII sources. As mentioned 
in Section 13.21 KDA approximates the density profile of the 
IGM with a modified King model (p = pQ(r/ao)~^). Here, we 
examine three density profiles and their affect on the maxi- 
mum lifetime of FRII sources. Two profiles (those of X[ and 
X2) are the parameters used in the KDA and BRW models for 
radio source evolution. We also implement an em pirically de- 
rived density profile of the IGM in local groups (iJetha et aTl 
I2OO7I X3 in Table|3]l. 

Our first density profile is associated with the default pa- 
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rameters of the KDA model. The average density po at the 
scale radius aq is 7.2 x 10"^^ kg m~^, o r less than 1 hydrogen 
atom per cubic meter ICanizares et alJ (|T987) argue that these 
values are typical out to a radius of 100 kpc from the source 
core. The value of /3 = 1.9 is larger than that of the other 
density profil e s disc ussed here but is still within the expected 
range. iFalld (Il991h showed that for /? > 2 the jet will not 
form shocks and thus no FRll sources will be present. KDA 
use this choice of beta as it agr ees with early X-ra y observa- 
tions of galaxies at low redshift (ICotter et al.lll996i) . We adopt 
po = 7.2 X 10"^^ kg cm"^, ao = 2.0 kpc, and /3 = 1.9 as our 
fiducial density profile. As noted previously, we calculate a 
Ufetime of 1 .22(±|j^) x 10^ yrs for Xis and 1 .15{±ll]) x 10^ 
yrs for Xib (Figure|5] left panel). 

In addition to KDA, the BRW model has also been used 
with success to explain the evolution of FRJl sources with 
time. The BRW model adopts values of po = 1 .67 x 10"^^ kg 
m"^ ao = 10 kpc, and /3 = 1 .5. iGarrington & Conway! (1 19911) 
use the depolarization of polarized synchrotron emission from 
the lobes of radio galaxies to suggest that radio galaxies pref- 
erentially reside in poor group environments. BRW ascertain 
a density profile from ROSAT ob servations of such groups 
jMulchaev & Zablud off 1998; Will ott et al.|[T998l) . Again, the 
average density at the scale length is less than one atom of 
hydrogen per cubic meter It stands to reason that this aver- 
aged, measured density profile is a good representation of the 
IGM in poor groups. We therefore perform our MLE using 
the BRW density profile (A2). Our likelihood as a function of 
tmoj: is shown in the middle panel of Figure|5] The most likely 
t„,, is 1.15(±»:?2) X 10^ yrs {Qs, X2s) and 1.59(±»:i^) x 10^ 
yrs {Qb, '^2b)- These values are in good agreement with those 
obtained with the KDA density profile. To good approxima- 
tion, the KDA model depends upon the quantity poAq rather 
than independent variations in each term. The BRW density 
profile product is a factor of five less than that of KDA yet 
only changes imax by a factor of ~ 0.8. 

With the advent of high resolution X-ray observations with 



Chandra, gas densities in galaxies and groups are being pre- 
cisely measured as a function of distance from the galaxy 
core. While one would like to know the distribution of gas 
density profiles seen in groups and clusters, the best avail- 
able data a re de nsity measurements of group galaxies by 
I Jetha et aTl (l2007l) . Using Chandra, they looked at 15 nearby 
groups of galaxies and computed a best fit power law to the ra- 
dial gas temperature and density profiles of these groups. We 
compute the maximum lifetime for this density profile (pa- 
rameter set A3 in Table Is) and find imax= 8.8(±og) x 10^ yrs 
for Qs and 1.31(±J];[^) x 10^ yrs for Qb (see also Figure |5l 
right panel). The Jetha density profile product is 24 times 
smaller than our default and alters t„,af by a factor of ^ 0.7. 
The similar maximum ages calculated with these three den- 
sity profiles strongly suggests our results are a weak function 
of IGM density profile. We further explore this behavior in 
the Appendix. 

4.2.2. Axial Ratio 

The axial ratio characterizes the shape of the lobe. More 
collimated lobes reach longer distances in a shorter amount of 
time as the area of their hotspot is smaller, allowing the jet's 
ram pressure to build up faster. We note that \ .Q <Rt < 6.0 
have been found for FRll s ources jL eahv & Williams 19841 
Leahy et al.lll989 ). although lMachalski et all (2004) gathered 
a sample of radio sources from the literature and found axial 
ratios as high as 7.5. However, this sample was biased to- 
wards the most powerful radio sources known and may not 
be indicative of typical FRJls. iKaiser & Alexander (1999j) 
found that Rj is dependent upon j3 while BRW and others 
postulate that axial ratio is a function of jet power Either sce- 
nario would help expla in some of the high axial ratios found in 
iMachalski et al.l (12004 ). As the vast majority of FRJl sources 
appear to have axial ratios between 1.0 and 6.0, we focus our 
attention on this region of parameter space. 

In addition to our default value of Rt= 2.0, we also inves- 
tigate Rt= 1,4, and 6. Our C{\.max) for all eight of these pa- 
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Fig. 6. — Maximum age t,„ax as a function of Rr (solid symbols and lines, 
left ordinate) and p (open symbols and dashed lines, right ordinate). Except 
for those pai'ameters that ai'e explicitly varied to test their degeneracy with 
t,„av, default values of the KDA model input parameters ai'e used. Both the 
Qs (squares) and Qg (triangles) jet power distributions are shown. Note that 
tmn.v shows no dependence on p. Rj has a strong effect on t,„a.x but is still 
within a factor of ~ 2 of the default value even at the extremes of the tested 
range. 

rameter choices is shown in Figure |6] When Rj decreases, 
the pressure at the hotspot increases as the lobe encounters a 
higher resistance from the IGM due to the increased surface 
area of the cocoon (KDA). Consequently, the hotspot advance 
speed decreases and we therefore expect the FRll lifetime to 
be inversely related to Rt- If the lobes are less coUimated, 
sources must be older to reach a given linear size. We refer 
the reader to Table[3]for our computed values of the FRll life- 
time for T??- = 10,2.0,4.0, and 6.0. 

4.2.3. Injection Index 

The injection index p (Section 13.2b is used solely in the 
calculation of the source luminosity. We refer the reader to 
KDA for an explanation of the relationship between luminos- 
ity and p. Here, we test the stability of our results against 
changes in the power law describing the injection of energetic 
particles into the jet. We have noted that substantial theoret- 
ical and observational evidence exists for 2.0 <p< 3.0 (Sec- 
tion l3.2b . To sample this range, we examine p = 2.5 and 2.9 
in addition to our fiducial value of p = 2.14. Recall that the 
lifetime for FRll sources under our default conditions {p = 
2.14) was 1.22(±°;2^) x 10^ yrs (Xis) and in5{±l^) x 10^ 
yrs {X\b)- Upon changing the injection index to 2.5 {X-is), 
the most Ukely t„„„- is 1.24(±*5;^^) x 10^ yrs; p = 2.9 (Xf^s) 
yields a lifetime of 1.26(±y jg) x 10^ yrs (Figure |6]l. Using 
the Blundell jet power distribution, we measure lifetimes of 
1.79(±^:i^) X 10^ yrs {X^B) and 1.79(±!j:Jf) x 10^ yrs (A-g^). 
While variation of p should not greatly influence radio lumi- 
nosity at a given frequency, we have confirmed that changing 
the value of p within the theoretically accepted range pro- 
duces no significant change in the lifetime estimate. 

We assert that our results are pseudo independent of source 
luminosity due to FIRST'S sensitivity at this observing fre- 
quency and the expected luminosity range of FRJls. Unless 
other models compute source luminosity in drastically differ- 



ent ways, we expect our results to be robust to changes in 
radio source evolution model. In particular, as the models of 
BRW and MK are based upon the KA model, use of either of 
these models should produce no appreciable changes to these 
results. 

4.2.4. Summary 

We have measured the range of t^ax consistent with our ob- 
served sample given the uncertainties in the parameters that 
characterize the parent population. Estimates of t^ax vary 
from 5.3 x 10* to 3.55 x 10^ yrs. Within the accepted range 
of input parameters described above, the most uncertain vari- 
able is Rt- Extreme lobe axial ratios, compared to our fiducial 
value, produce a factor of 2 change in t^ax- 

While one could allow for a distribution in density pro- 
file and axial ratio as part of this analysis, we judge that 
our present sample size is too small to provide meaningful 
constraints. Instead, we have examined a wide and reason- 
able range of both density profile and axial ratio. If distribu- 
tions (as yet unquantified by observations) in these parame- 
ters were used instead, the resulting t„iax would be a weighted 
mean of the results already obtained here. When generally 
discussing FRll lifetime in subsequent sections, we conserva- 
tively choose t„a^= 1.5(±0.5) x 10^ yrs. Other choices of A", 
do not alter the lifetime by more than 20% compared to those 
found with Xis and Xib- 

4.3. FRIIDuty Cycle 

In addition to their lifetime, the duty cycle as sociated with 
FRJls is an important timescale. Using the B erlind et al.l 
( ,2006 a) catalog of groups, we measure the typical time be- 
tween episodes of FRll radio activity in central group galax- 
ies. Of the 2020 groups studied, 26 central galaxies were host 
to FRll radio sources. We determine the level of incomplete- 
ness from our mock catalogs. The default simulation with Qs 
{X\s) yielded a lifetime of 1 .2 x 10^ years. Of the 10* sources 
in X\s, 6x10^ would have been detected by FIRST and iden- 
tified as FRJls. We therefore estimate that our observed sam- 
ple of FRll sources is 60% complete for sources that follow 
the quoted model parameters (e.g. those with initial jet pow- 
ers above Qmin)- Sources may not be detected for two reasons: 
(1) Too young and therefore smaller than the angular resolu- 
tion of FIRST or our required angular size to be considered an 
FRll and (2) Luminosity corresponds to a flux below FIRST'S 
sensitivity at the source redshift. We conclude that the incom- 
pleteness corrected FRll duty cycle is 2.2% (44/2020) in cen- 
tral group galaxies. The corresponding average time between 
FRJl phases in these central galaxies is 5.6 x 10** years. Due 
to a different completeness level for Qb {X\b), 1 .7% of central 
galaxies exhibit FRJl activity using this jet power distribution. 
The time between FRJl phases is thus a factor of 2 higher at 
1.0 X 10^ years. We note that there are many FRJ sources in 
this group sample as well and consequently the environments 
of some of these sources may not be amenable to an FRll mor- 
phology. This implies the FRJl duty cycle is actually higher 
for those group environments suitable for FRlls. 

5. RESULTS: POPULATION STATISTICS AND CLUSTERING 

Our attention was initiall y brought to the Berlind group 
sample dBerlind et al.l2006al) by their subsequent paper on the 
clustering or bias of grou ps as a function of v arious group and 
central galaxy properties dBerlind et al.l l2006b). It is interest- 
ing to determine if these properties are correlated in any way 
with radio activity, and whether radio activity is correlated 
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with m ore clustered systems. Here we quickly outline the re- 
sults of iBerlind et alJ (l2006bl) and discuss how radio activity 
varies with a number of the parameters studied therein. 

Group an d central galaxy prope rties tested for correlations 
with bias in iBerlind et al. ' (2006b) were group number rich- 
ness Ngrp, group velocity dispersion a,,, total group color (g- 
T),ot, central galaxy color (g-r)^™, central galaxy luminosity 
M, c(.„, and the magnitude gap between the central galaxy 
and the second brightest member of the group AMr- As- 
suming group luminosity and mass were monotonically re- 
lated and matching t he group LF to the ha lo mass function of 
IWarren et all (l2006h . IBerlind et aP (l2006bh divided their sam- 
ple into four bins of group mass centered at 10'^^, 10'^ *\ 
lO'^'^, and lO''**' h"'M0. For every galaxy or group pa- 
rameter, each mass bin was split into "high" and "low" sub- 
sets, based upon the median parameter value in a given mass 
bi n. The ratio of the hig h and low bias functions (see eq. 3 
of iBerUnd et al.l ('2006b')) was then evaluated for each prop- 
erty across the mass bins. They found that only Ngrp and ct,, 
were strongly correlated with clustering over their entire mass 
range and conclude that bluer galaxies are located in more 
strongly clustered groups than redder central galaxies. For 
the most massive systems, they also found AM,, is a function 
of bias while they reported no significant correlation between 
clustering and M^cen- In light of these results, we examine the 
trends between radio activity and each of these parameters. 

Our sample of 2020 central group galaxies is comprised of 
the most luminous groups in the Berlind catalog. We deter- 
mine the radio fraction of central galaxies (r/rac) with three 
distinct definitions of what constitutes a radio source: the cen- 
tral galaxy simply contains a core source or is host to extended 
radio emission (either FRI or FRII) (r/raci)', the radio source 
has LpiRST > 10^2 5 W Hz-1 (rfraci) as in lBest et al.1 ( l2006l) : 
and most conservatively, the cent ral galaxy is radio loud with 
logFfjRsr/iogFsDss > 1 (rfraci, llvezic et al.ll2002h . By any 
of these definitions, a larger rfrac corresponds to a greater like- 
lihood of harboring an active radio source. We broadly expect 
larger, more clustered systems to have a larger viral radius and 
possibly more virialized gas than smaller systems. If radio 
mode feedback is indeed a mechanism to combat the cooling 
flow problem, one might expect a correlation between rjrac 
and bias. 

We first determine the radio fraction as a function of group 
richness. Figure [T] (upper left) shows r/mc as function of 
logNgrp- Roughly 20% of the central galaxies are host to ra- 
dio sources even in the smallest systems and there is strong 
evidence for a rise in rjrac with increasing Ngrp- While this 
behavior is only at the one to two a level for rfrac,i, the slope 
of the rjrac-^ogNgrp relation increases for rfrac,i and rfrac,3- 
At the 3.5(7 level for rfrad, groups of 10 or more member 
galaxies are almost three times as likely to harbor a radio 
loud source in their central galaxy than smaller groups. Even 
though this trend may suggest even higher rfr„r in br ight- 
est cluster galaxies (BCGs), work by iBest et alf ( l2006l) de- 
termines a rfrac maximum at approximately the value we find 
(35%). 

Two quantities used as an indicator of group mass are ct,, 
and Mrjot- While we are hindered by low number statistics 
in the highest bins for ct,,, we see no statistically significant 
trend of r/mc with the velocity dispersion of the group. There 
is marginal evidence for a correlation between r/rac and M, /o,, 
though not enough to be conclusive. 

The central galaxy's properties may be more important than 



the group's for the radio fraction. To investigate, we look 
for any correlation between r/rac and Mr.cen, (g-r)ccn, and 
AM,- (Figure IT). The luminosity of the central galaxy is the 
strongest indicator (at the 5(7 level) of radio activity that we 
studied (Figure |7] upper right). When compared with the 
faintest central galaxies in our sample, r/rac is increased in 
the most luminous galaxies by a fact or of ~ 12. This rel ation- 
ship had been shown previously by iBest et akl (l2005a) who 
found that the radio-loud fraction was a strong function of 
host galaxy mass {fradw-ioud oc M^'^). If we assume that the 
mass to light ratios of the brightest group galaxies do not vary 
wildly, this work demonstrates the continuation of this trend 
towards lower mass central galaxies. 

The color of the central galaxy {g-i)cen may also indicate 
the presence of a central radio source as this parameter does 
correlate with clustering properties. We find an increase at 
the ^ 3(7 level from ~ 20% to ^ 35% in rfrac,\ in the red- 
dest central galaxies (Figure]?] lower left). The effect is even 
more pronounced when looking at the radio loud fraction via 
'"/rac.2 and r/rac.i and therefore redder, presumably older ellip- 
tical galaxies, have an increased probability of hosting a radio 
source. We do not find a significant AGN contribution to the 
spectra of our FRII sample and do not expect any systematic 
color enhancement due to AGN. Most of the central galax- 
ies are quite red, yet it is interesting that the reddest galaxies 
are most prone to radio activity. Berlind et al. (2006b) found 
that bluer galaxies were in more clustered environments, es- 
pecially in the high mass regime. Their results combined with 
our own suggest that radio activity and clustering bias are not 
necessarily correlated. 

Central galaxy identification with FoF gro up and clus- 
ter find ing algorithms is not perfe ctly robust. [Berlind et al.l 
(l2006ab and the maxBCG method (iKoester et al.ll2007h bofli 
rely upon the observed correlation between galaxy luminos- 
ity and radial distance from the center of the potential well 
to designate central galaxies. The scatter in the luminosity- 
position relationship decreases with increasing luminosity gap 
between the brightest and next brightest group or cluster 
member (AM^). We find a positive trend between AM^ and 
the radio fraction (Figure |7] bottom right). As AM^ increases 
from to 1.75, radio fraction increases ^ 3.5, 7, and 4 fold 
iffrac,\, ffmc,2, ffrac,3, respectively). The luminosity domi- 
nance of a galaxy in the group environment, and hence its 
certainty as the central galaxy, is strongly correlated with ra- 
dio emission. 

6. SUMMARY AND CONCLUSIONS 

We present the first estimate of the maximum radio source 
lifetime in a well defined sample of groups and clusters 
of galaxies selected from the SDSS and FI RST surveys. 
Cross-correlating the group catalog identified by lBerlind et al.l 
(2006a) with FIRST, we identify FRII sources associated with 
central group galaxies. The projected lengths of these sources 
are determined by the separation of the end of the lobe and 
either a core radio source or the optical position of the host 
galaxy. From this sample and measurement, we determined 
the most likely maximum lifetime of FRJIs through compar- 
ison with mock catalogs of known jet, lobe, and IGM prop- 
erties. We modeled the linear growth of these mock sources 
with time and their luminosity evolution with the prescription 
set forth in KA and KDA. By fixing the jet, lobe, and envi- 
ronment properties, the maximum amount of time the jet is 
fed energetically by the AGN was determined via a maximum 
likelihood approach. Our prototype for the average state of 
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FRII lobes and surrounding medium is listed in Table[3](-^i). 
We stress that our sample is representative of the group en- 
vironment and as such our lifetime estimate is for a typical 
FRJl source in these systems. While individual sources of 
lengths and lifetimes exceeding those seen here are observed, 
they are far less numerous than the members of this sample. 
Using the default input parameters to the KDA model and 
the Blundell Q distribution, our sample's most likely t„Mx is 
\.15(±l}^) X 10^ yeai-s. 

We numerically investigated the dependence of any single 
input parameter and t,„ai, while we present an analytic treat- 
ment in the Appendix. First, the power associated with the jets 
can affect the lifetime calculated in our analysis. We use two 
distributions: (1) One that repro duces the act i ve rad io galaxy 
luminosity function derived in ISadler et"aD (l2002h and (2) 
The FRII jet power distribution found on energetic grounds 
from BRW. Keeping all other input variables fixed, we find 
that changing the jet power distribution from Qs to Qb pro- 
duces only a ^ 30% change in the value of t,„,„ (1 .22 x 10^ to 
1.75 X 10^ yr). Mock radio sources were expanded into three 
distinct IGMs. In addition to the KDA density profile, the less 



dense profile of the BRW model a nd a recent observational 
measurement by Jeth a et all (l2007h were tested. As expected 
the calculated lifetime decreased relative to our fiducial t,„ai, 
but only marginally so: ~ 10% for BRW and ^ 30% for the 
Jetha profile. In practice, the density profile of the IGM in 
galaxy groups spans some distribution whose shape will con- 
tribute to a range of expansion speeds for FRII sources. The 
work here demonstrates that t„„u is a weak function of envi- 
ronment and will likely vary by less than 50%. We find that 
the axial ratio of the lobe had the largest effect on the cal- 
culated lifetime of our observed sample. Still, the calculated 
lifetime of extreme axial ratio lobes was not different from 
our fiducial results by more than a factor of ^ 2. We sam- 
ple the entire range of expected axial ratio and find t„„u in the 
range of ~ 0.5-2.4 x 10'' years {Qs)- Due to the sensitiv- 
ity of FIRST and the relatively small distance to the groups 
in this study, uncertainties in the particle energy distribution 
do not effect our results. We created mock catalogs that only 
differed in the choice of the injection index p and found t„„v 
was statistically identical over the entire probable range. The 
values and associated uncertainties of parameters solely con- 
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tributing to the luminosity of the sources are insignificant in 
this analysis. 

Taken together, these tests on the parametric dependence of 
t,„fl v demonstrate how robust our results are to ambiguity in the 
true distributions of many jet, lobe, and environmental param- 
eters. Due to the self-similar nature of the length to age rela- 
tionship determined in KA, these results would not change ap- 
preciably if another radio source evolution model was incor- 
porated into our analysis. While individual extreme sources 
may be observed with ages significantly larger than the maxi- 
mum lifetime found here, the vast majority of FRll sources re- 
side within the parameter space examined. We therefore con- 
clude that the average FRIl lifetime is 1.5(±0.5) x 10^ years. 

This characteristic maximum lifetime of the FRll popu- 
lation is in good agreement with various methods of radio 
source dating for individual objects. Recently, O'Dea (2007, 
in prep.) used spectral aging to estimate the source ages of 31 
powerful FRll radio galaxies. Under minimum energy con- 
ditions, all source ages were less than 1. 1 x 10^ years while 
most were several Myr old. Allen et aLl (12006^ presented a 
comprehensive study of luminous elliptical galaxies, measur- 
ing the age of (primarily FRJ) radio sources still energetically 
associated with their parent AGN. They estimate the cavity 
age using the sound speed of the IGM and the size of the bub- 
bles in the X-ray gas and find ages of 10'' to 10** years. 

The maximum lifetime of 1.5 x 10^ years for FRJls is also 
in good agreement with observational constraints on the life- 
time of_QSOSj_which suggest values from 10^ to 10** years 
(e.g. lMartinill2004l) and the value of the e-folding or Salpeter 
timescale for black hole growth at the Eddington rate = 
4.5 X 10^ yr (for 10% radiative efficiency). This is most likely 
a coincidence, because these FRlls show little evidence of nu- 
clear activity in their SDSS spectra, are consequently accret- 
ing at a low rate, and also show no evidence that they have 
been triggered by mergers. In fact the detailed models of 
AGN triggered in major mergers find that the lifetime is pro- 
gressively longer when the AGN is defi ned at a progressivel y 
lower luminosity or accretion rate limit (Hopkins et al.ll2065h . 
At the accretion rates of these FRlls, the predicted Ufetime 
would be substantially longer than 10^ years in the merger 
scenario. 

We investigated the connection between host galaxy and 
group properties and the probability of harboring a radio 
source. The size of a group correlates with r/„c at the 2a 
level. Assuming some average galactic mass, group size 
should be proportional to its viral radius. It is plausible that 
larger viral radii would encapsulate greater gas mass on aver- 
age. If AGN feedback regulates star formation, larger systems 
would need shorter duty cycles or longer lifetimes of feed- 
back activity to counteract their stronger cooling flows. The 
increase in r /mc with group size supports the idea of episodic 
AGN heating. 

The luminosity of the host galaxy is the strongest predic- 
tor of radio activity. Comparing galaxies with Mr = -20.5... - 
23.0, the radio fraction increased with luminosity by at least 
a factor of 5 at the greater than 4a level depending on the 
criterion used to determine the radio source activity thresh- 
old. This result is in agree ment with previous studies of radio- 
loud AGN in BCGs (Bes t et al.ll20"06HCroft et al.ll2007l) . Our 
combined results point toward galaxy luminosity being the 
most important indicator of radio activity over a large range 
group/cluster mass. A slightly weaker correlation (~ 3a) 
between host galaxy color and rfmc was also found. Red- 
der galaxies have an increased probability of harboring ra- 



dio sources. Interestingly, galaxy color was found to be 
anti-correlated with c lustering bias recently in the literature 
jBerlind et al.ll20 06F). The two results suggest no connection 
between radio fraction and bias. 

The luminosity gap between the brightest and second 
brightest group or cluster member has been discussed as a 
measure of certainty regarding c entral galaxy identific ation. 
The maxBCG catalog and that of ' BerUnd et al.1 (l2006as) both 
designate the brightest cluster (group) member as the central 
galaxy. A caveat of this technique is uncertainty in the cen- 
tral galaxy's position with respect to the center of the group 
or cluster's potential well. As A M^ increases, the po sition- 
luminosity relationship tightens (iLoh & Straussll2006l) . We 
find a significant correlation between r/rac and AM^, which 
may be due in part to misdiagnosis of the central galaxy or 
that this g alaxy is not at the center of the group potential. 
iBest et alT (12006) examined the radio loud AGN fraction of 
BCG galaxies, finding an increased likelihood of radio source 
association when compared to field galaxies. The discrep- 
ancy in radio fraction between central and field galaxies may 
increase if the possible misclassification of central galaxies is 
taken into account. This potential source of incompleteness 
associated with all central galaxy studies should be addressed 
in future analysis. 

Our results have important implications for AGN feed- 
back models. The recent literature includes many mod- 
els for AGN energy injection into the ICM and its sub- 
sequent effects on gala xies and hot, gaseous halos (e.g. 
ISijacki & Springell2006l) . Heating from AGN is now included 
in many structure formation models to match the observed 
galactic LF (e.g. Crotonetal. 2006; Springel et al. 2005a,lJ 
iHopkins et al.l l2006l) . The FRJI lifetime presented here is 
an accurately measured timescale over which jets are ener- 
getically sustained by AGN and perform work on the ICM. 
Many models use an instantaneous injection of energy into 
the ICM as the timescales involved have historically had or- 
der of magnitude uncertainty. While it is true that FRIs are 
far more numerous than FRIIs, there is increasing evidence 
of FRIIs and FRIs being the manifestation of the same under- 
lying jet mechanism in different environments. A dramatic 
empirical demonstration of this are the hybrid morphology 
sources (HYMORS) in which one lobe of a double source ap- 
pears to be an FRII while the other exhibits FRI charact eristics 
(iGopal-Krishna & Wiitall2()00l: iGawronski et al.ll2006h . If FR 
morphology is a function of environment, it is reasonable to 
assume both classes share the same lifetime. We have shown 
the lifetime of FRJI sources to be ~ 10^ years in a wide range 
of initial conditions and it is most accurate to inject energy 
over this timescale. 
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APPENDIX 
PARAMETER DEPENDENCE 



Here we derive the analytic dependence of t„,aj on the other KDA model parameters. By plugging eq.[3]into eq.|2]we obtain 
the full equation for the length of a lobe in terms of the parameters discussed in section [ 



I = cip,RT) {pa^y'"'~^^ gl/(5-/3)f3/(5-ffi^ (Al) 

where c{(3,Rt) is just ci from eq. |2]and eq. 4 (KDA), but shows that it is dependent upon /? and Rt (see eq. 25 of KDA for 
the definition). We can easily manipulate this equation to solve for the age of the source, t. To zeroth order, we are trying to 
determine t„„„ such that a source at the average age f = has the median lobe length I. Hence setting t = t and I = I we find 

t = p-f'^/'cipr-'^' (pa^) Q-"\ (A2) 

Note that t and therefore f^ax only weakly depends on the density distribution, pqQq, and is independent of p. 

We first compare two different density distributions with pi, a^, /3i and p2, 02, Pi, respectively. The distribution of jet powers 
is kept constant (either Sadler or BRW). In this case the two maximum ages are related by 



fmax.2 



f-j r(5-/32)/3 c(/?2)(ft-5)/3 



(A3) 



Assigning the KDA density profile (^"1) to tmax,u the density profile from iJetha et al.l (l2007l) {X{) to f,„nA.2; and letting 
h = 25.75 kpc {Xys) and h = 30.85 kpc {Xis), we find t„ax,KDA.Qs I t„,axje,im,Qs = 1-36 and with a similar analysis tmax,KDA,QB I 
tmaxjetha,QB = 1-35. Thcsc rcsults are virtually identical to those found in our statistical approach, namely tmax,Xishmax,Xzs - 1-39 
and t,nax,XxBhmax,x^B = 1-34. Similarly, analytically comparing the KDA and BRW density profiles, t,„ax.KDA,Qs I tmax,BRW,Qs = 1-06 

and tmax,KDA,QB I tmax,BRW,QB = l-^O whilc tmax,Xis / tmax,X2s = 1-06 and t„,ax,XiB / tmax.XjB = I-IO- 

To isolate the effect on t^ax by the choice of jet power distribution, we hold the density profile and axial ratio of the lobes fixed. 
In this case, the ratio of t^ax simply becomes a function of the Q distributions: 

- ,-(5-/3,)/3 . . -1/3 

'^max.sad _ had _ '1 { \C-sad \ /'A/l\ 



Taking the median value of each Q distribution as a representative jet power and J from our mock catalogs, we find 
tmax.KDA.Qs / tmax.KDA,QB = 0.69. Again, this analytical result is in excellent agreement with the simulation derived t„ax.Xis/fmax.XiB = 
0.70. Notice that the statistically obtained ratio of Qs to Qb across all parameter choices in tableware between 0.67 and 0.72. A 
similar analysis for Xis / X/b reveals that the theoretically and practically obtained ratios are in very good agreement. An increase 
in the median of the Q distribution by a factor of 2.15 (from Qb to Qs) decreases tma, by only a factor of ^ 0.70. 

Finally, we examine the relationship between axial ratio and t„a^. As ci and I are the only terms in equation |A2I that depend 
upon Rt, the ratio of lifetimes derived from two models that only differ in their choice of Rj is: 

wx.«.., r/'-'''^/^(/3i)(ft-5V3 

(Aj) 



Wx,«.. f«„ z-(5-A)/3c(/32)(A-5)/3- 

As an example, w e fin d t,nax(RT= 2.0,2s) / t,„„v(/?r= (!>.Q,Qs) = 2.33. This expectation is met very well as the ratio of Xis I Xes 
= 2.30. In fact, eq. lA5l is a good fit to the \.,„ax distribution resulting from simulations that vary in Rj- 

The equations here derive only from the self similar solution for the length of these FRll sources (eq. |2|. This particular 
equation is a characteristic of all the radio source evolution models discussed in this paper (KA, KDA, BRW, MK). We have 
shown that \.,„ax is a weakly varying function of the density profile, jet power distribution, and axial ratio. The true nature of 
these jets, lobes, and their environments must be systematically and significantly different from the assumptions adopted here to 
change our results by a factor of 3 or higher. 
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